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Abstract. We present a differential measurement of the azimuthal anisotropy of 

charged hadron production in Pb+Pb collisions at v /snn=2.76 TeV. This azimuthal 

anisotropy is expanded into a Fourier series in azimuthal angle, where the coefficient 

ij_) ' for each term, v n , characterizes the magnitude of the anisotropy at a particular angular 

scale. We extract v 2 — Vq via a discrete Fourier analysis of the two-particle A0 — A77 
correlation with a large Arj gap (| Arj\ > 2), and via an event plane method based on the 
Forward Calorimeter. Significant v 2 — vg values are observed over a broad range in px, 
r\ and centrality, and they are found to be consistent between the two methods in the 
transverse momentum region px < 3 — 4 GeV. This suggests that the measured v 2 — v% 
obtained from two-particle correlations at low px with a large A77 gap are consistent 
with the collective response of the system to the initial state geometry fluctuations, 
^1 \ and is not the result of jet fragmentation or resonance decay. 



CN 



o 



Among the many striking results obtained at RHIC, one important observation is the 
novel "ridge" -like and "cone" -like structures of the two-particle correlation in relative 
azimuthal angle A0 = 4> a ~4>b an d pseudorapidity A77 = r] a — r] h in Au+Au collisions p] . 
These structures, obsent in elementary proton-proton collisions, are found to extend 
over a large Ai] range, and show rich p^ and centrality dependence. They were initially 
interpreted as response of the medium to the energy deposited by the quenched jets [2]. 
However, recent studies j3] suggest that they could be related to the initial geometric 
fluctuations and strong collective flow. In this scenario, the spatial fluctuations of 
nucleons lead to shape deformation at various angular scales, which induce high-order 
anisotropies of the emitted particles through collective expansion: 

00 
dN/d<p oc 1 + 2 Y, Mpt) cos (n (<p - tf n )) , (1) 

n=l 

where v n is the magnitude of the n th order harmonic flow [4J. The two-particle 
correlation, being a simple convolution of two single particle distribution, is naturally 

influenced by the same harmonic flow: 

00 
dN/dA(p oc 1 + 2 J2 v n (p T a )v n (p T h ) cos (nA<f>) (2) 

n=l 
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where the phase of n th harmonic flow \l/ n (known as the event plane or EP) drops out 
in the convolution. 

The v n can be extracted directly by correlating single particles with the measured 
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order EP (the event plane method) or via a Fourier transform of the two-particle 
azimuthal correlation (the two-particle correlation method). We present v n results from 
both methods based on ~ 8 yub -1 Pb+Pb data from the 2011 LHC heavy ion run at 
v /snn = 2.76 TeV [5]. A detailed comparison is made between the two methods, followed 
by an interpretation of the long-range structures in terms of geometric fluctuations and 
collective flow. 



1. Event Plane Method 

The event plane \l/ n is estimated using the transverse energy flow measured in towers 
from the Forward Calorimeter (FCal) within 3.3 < \rj\ < 4.8. By default, each v n is 
measured by correlating tracks in the Inner Detector covering \rj\ < 2.5 with the EP 
from the full FCal, followed by a resolution correction [1]. In a slightly modifed method 
(FCalp(N)), tracks are correlated with the EP calculated with FCal in the opposite 
hemisphere, i.e tracks with r\ > (rj < 0) are correlated with the EP obtained from 
FCalN (FCalp). The FCalp(N) method greatly increases the pseudorapidity separation 
between the tracks and the EP from about 3 units to about 5 units, thus it is 
much less affected by the so-called "non-flow" correlations, which stem primarily from 
jet fragmentation and resonance decays. This is especially useful for measuring the 
pseudorapidity dependence of v n . We present results from both methods, relying more 
on the full FCal method when better resolution is needed (e.g. for the higher-order 
harmonics) . 

Figure Q] shows the resolution factor in 5% centrality intervals along with a 0-1% 
most-central interval from the full FCal method. Significant values are observed for n=2- 
6. It is interesting to note that the resolution factor for n = 3 exceeds that for n = 2 in 
the 0-1% most central collisions. Figure [2] summarizes the r\ and pp dependence of v 2 — v e 
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Figure 1. The resolution (left) and resolution parameter Xn 14! (right) vs centrality 
(smaller number refer to more central collisions) for n — 2 — 6 from the full FCal 
detector, together with their systematic uncertainty bands. 
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Figure 2. (left) r\ dependence of v n for 2 < px < 3 GeV from FCalp(N); (middle) px 
dependence of v n ; (right) px dependence of v n /v 2 , for two centrality selections. 



for two centrality selections. A very weak rj dependence (less than 5%) is observed for 
all n within \q\ < 2.5. 

All measured harmonics have a similar px dependence: they first increase rapidly 
up to pt ~ 3 — 4 GeV and then fall. The similar px dependence motivates us to 
find a simple scaling relation between different v n : v\[ n oc v 2 ■ This is qualitatively 
understood in a "blast wave" scenario, where v n oc (3 n , with /3 being the radial flow 
velocity ^j. However, it should be pointed out that this scaling is only approximate. In 
particular, the ratios for n = 4 — 6 are close to each other, while they are systematically 
higher than those for n — 3. 

In Fig. [2l v 2 — v 5 are shown only up to 12 GeV, so the low px region can be 
seen more clearly. However, ATLAS has measured v-i out to much larger px, as shown 
in Figure |3] [7], with much smaller statistical errors at high px than previous RHIC 
results [8]. The v 2 clearly continue to drop out to 10-12 GeV, and only vary slowly 
beyond that, but at a level that is consistent with pQCD calculations [9]. 



2. Two-particle Correlation Method 

The correlation function is constructed by dividing the same-event pairs with mixed- 
event pair, with a pair acceptance extending to |Ar/| = 5. The normalization is fixed by 
equating the counts of the same-event and mixed-event pairs in 2 < | Arj\ < 5, which is 
then applied for all Ar] slices. Each 1-D correlation function (obtained by integrating a 
selected | Arj\ range) is expanded into a Fourier series, with the coefficient t> n , n calculated 
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Figure 3. v 2 vs pt at mid-rapidity for several centrality selections. 



directly via a discrete Fourier transformation (DFT): v n>n = (cos(nA0)). Figure @Ji 
shows one such projection for 2 < \Arj\ < 5 and the corresponding contributions from 
individual harmonic components. 

If the observed modulations are due only to collective flow, then we expect v n>n to 
be factorizable into the product of two single-particle flow coefficients: 

V n ,n(PT^PT h ) = MPT^MP^) (3) 

Thus for correlations where the two particles are selected from the same pp range, we 
calculate the single-particle harmonic coefficient as v n = y/v n ,n- One such example 
is shown in Figure Hb. We have repeated this procedure for each At] slice and the 
results are shown in Figure Hb. The peak at small At], which comes from near-side jet 
fragmentation and resonance decay, is excluded by requiring 2 < \Arj\ < 5. 



■e- 
< 



ATLAS Preliminary 
2<|AT!|<5 
2<p", p£<3 GeV 




1(T 



0.005 



-0.005 



b) 



ATLAS Preliminary 

J Ldt = 8(ib~ 1 

2<|Ati|<5 



iiLl 



i IM i i i 1 i i i i 

-e- V. 

-B-V] 



> 

0.05 



t'tHMj 



'I 1111 ! 1111 ! 1111 ! 1111 ! 1 

ATLAS Preliminary 
I Ldt = 8fib' 1 
2 < p",p£ < 3 GeV 
0-5% 



"S^ooo o o o o o 
•••• • • * • 

fan nn nan a □ n D n 



=*v 



o o o o o o 



*** o * o a * 



A<|> 



AT| 



Figure 4. The steps involved in the extraction of the v n (2-3 GeV fixed-px correlation 
in 0-5% centrality): a) AcA correlation function for 2 < |Ary| < 5, overlaid with 
contributions from individual Fourier components and the sum, b) Fourier coefficient 
v n vs n, and c) v 2 — v§ vs Arj. 



In order to confirm that the v n coefficients from the two-particle correlation method 
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Figure 5. v n (pT ) = """ "^i-'f/ ^ vs A?y for target px of 1.4 < px b < 1.6 GeV, which 
are extracted from four reference px a bins, supporting the factorization relation Eq. [3] 



reflect collective flow, the factorization relation in Eq. [3] has been explicitly checked with 
correlations of pairs of tracks from different px ranges. For correlations with | Ar/| > 2, 
this factorization indeed holds for t> 2 — v 6 at 5%-10% level for px < 3 — 4 GeV in the 
70% most central events, where away-side jet contributions are small. One example of 
such explicit check is shown in Figure Conversely, the factorization is found to break 
down for vi, which could be due to the fact that v\ from collective flow changes sign 
going from negative 77 to positive 77 [10]. However, similar effects can also be caused by 
momentum conservation effects, e.g. the recoil of the away-side jet in a di-jet system. 

Figures [6] shows the extracted Fourier spectra for two fixed-px correlations with a 
large Arj gap (2 < |At/| < 5) in 0-1%. v-i is the largest at low pt, but it becomes less 
than f 3 and V4 for tracks in 2 < px < 3 GeV. The rate with which the v n decrease 
with n were shown to provide important insight on the acoustic horizon and viscous 
effects HTl. 
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Figure 6. v n vs n from same-charge, opposite-charge and all pairs in 0-1% centrality 
bin for two fixed-px correlations: 0.5-1, 2-3 GeV. 



ATLAS Vn results 



.1 " " i " " i " " i ' i " i " " i " " i " " i " " i 
ATLAS Preliminary 1-2 GeV 



-|""|"" |llll|lll l|l 



TO 



rjTTTTJ- 



I 



Ldt = 8|xb" 1 Stat. Error only 



10" 1 



> 



B a a a b b ■ 

A A " - - „ „ 
V v 9 y 9 



U □ 



10" : 



-9- ^V 5 

I .... I .... I 



* £ * *r ^ 

+ * + « 



. I .... I .... I .... I . 



©A 



eV 



9 9 ® » » 9 



aH BBBBDn n 

a a A 

9 V v v v v v 



****£ 






-Open Symbols: correlation method 

forfixed-p and2<|Ar||<5 

Solid Symbols: full FCal EP method 

for |ti|<2.5 



1.2 



a. 1 
w 

0.8 



o o o 



onnopooQ 




^tf^- 



A 



■a- -a 



20 30 40 50 60 70 80 90 100 

1-Centrality [%] 



20 30 40 50 60 70 80 90 100 

1-Centrality [%] 



Figure 7. Comparison of v n obtained from the two-particle fixed-px correlation 
method (solid symbols) to the full FCal results (open symbols) for 1-2 GeV (left panels) 
and 2-3 GeV (right panels). The ratios between the two methods are shown in the 
bottom panels. 



3. Comparison of the Methods and Discussion 

Figure compares the centrality dependence of v n obtained from the two-particle 
correlation method and those from the EP method for two pp intervals. The agreements 
is within 5% for v 2 — v$ over a broad centrality range, but worsen to about 10% for v 5 and 
15% for vq. However, they are within the quoted systematic errors shown in previous 
sections. 

The consistency between the two methods implies that the structure of the two- 
particle correlation at large Arj can be largely accounted for by the collective flow. We 
check this explicitly by reconstructing the correlation function from v n measured by the 
EP method as: 

C(A0) = 6 2p (l + 2vll cos A + 2 E v n P '^n P,h cosnA(j)) (4) 

n=2 

where the b 2p and v-^ are the pedestal and first harmonic term from the two-particle 
correlation analysis, while the remaining terms are calculated from the v n measured 
using the EP method. Figure M shows excellent agreement between the measured and 
reconstructed correlation functions. This is more striking for the 2-3 and 3-4 GeV bin 
in the 0-1% most central collisions, where even a more pronounced "double hump" 
structure can be well reproduced. The V\ \ term not measured by the EP method, 
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also plays a significant role, but it is not sufficient to explain the near and away-side 
structures obtained in the correlation function. Figure also shows that the away-side 
"double-hump" is the result of the detailed interplay between the odd harmonics ( v 3 
and v 5 ) and even harmonics (^2,^4, and v e ). 
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Figure 8. Correlation function data compared with that reconstructed from v\ : \ from 
two-particle correlation and i>2 — vq measured from EP method in 0-1% ccntrality for 
different pt ranges. 



If the low pt correlation function with large At] gap is dominated by harmonic flow, 
an important question remains as to where in phase space the jet contributions can be 
observed. Figure M shows the centrality evolution of the 2-D correlation function for 
particles with 2 < p^ < 3 GeV. While central events show structures that are long range 
in Arj as well as a shorter-range jet correlation around At] = A<p = 0, moving towards 
more peripheral events show that these long-range structures eventually disappear, with 
clear jet-related peaks emerging on the away-side. Figure [10] shows the px evolution 
of Acf) correlations in the 0-10% most central collisions. A rapidity gap of |Ar/| > 2 
is required such that the near-side peak reflects mainly the "ridge" contribution. Its 
magnitude first increases with px to 4-5 GeV then decreases, reflecting the fact that all 
of the v n reach a maximum at 3-4 GeV as shown in Fig. [2j A narrow away-side peak 
emerges in the 6-8 GeV bin, which quickly dominates the correlation structure at higher 
Px- This away-side peak presumably comes from the fragmentation of the recoil jet. 
This figure qualitatively suggests that the transition in p^< from being flow-dominated 
to jet-dominated, for correlations of tracks from the same px range, happens somewhere 
around 6-8 GeV. Interestingly, this is the pi region where the single hadron suppression 
is the strongest [12] . 
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Figure 9. 2-D correlations for 2 < py a , pr b < 3 GeV as a function of centrality. The 
near-side jet peak is truncated to better reveal long range structures. 



4. Conclusion 

In summary, the higher-order harmonic coefficients v 2 — vq have been extracted both by 
correlating tracks with the event plane determined at forward rapidity and by using the 
two-particle correlation method with a large pseudorapidity gap (|Ar/| > 2). Significant 
v 2 — v 6 are observed and they are consistent between the two methods for px < 4 GeV. 
The t>2 — Vq are found to decrease only slightly with \rj\ in < \r)\ < 2.5. All v n exhibit 
similar p^ dependence, namely, they all increase with p^ to around 3-4 GeV and then 
drop for higher px- However, the higher-order harmonics show a stronger px variation, 
which is found to follow a simple scaling relation, v\[ n oc v 2 , varying only weakly 
with px- We find the main structures of the two-particle correlation at | Arj\ > 2 can be 
explained by v 2 — vq which largely reflect collective flow, and a v\ term which accounts for 
momentum conservation effects (possibly including a global directed flow). We conclude 
that the low px correlation functions do not allow significant contributions from a short 
Arj range medium response to jets. Fluctuations in the initial state geometry, along 
with a non-zero viscosity of the medium, are potentially responsible for the detailed 
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Figure 10. px evolution of two-particle A0 fixed-pT correlations for 0-10% centrality 
selection, with a large rapidity gap (| Ari\ > 2) to suppress the near-side jets and select 
only the long range components. 



behavior of the harmonic coefficients. A detailed comparison with viscous hydrodynamic 
calculations p3] can help elucidate the nature of these fluctuations and better constrain 
the transport properties of the hot, dense medium. 
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